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COMPUTER SIMULATION OF AIRCRAFT MOTIONS 
AND PROPULSION SYSTEM DYNAMICS FOR THE YF-12 
AIRCRAFT AT SUPERSONIC CRUISE CONDITIONS 
Stuart C. Brown 
Ames Research Center 


SUMMARY 


This report describes a computer simulation of the YF- 12 aircraft motions 
and propulsion system dynamics. The aircraft is a delta-wing twin-engine type 
which is capable of cruising at high altitudes and supersonic Mach numbers and 
the simulation is intended to represent the aircraft and its systems for these 
conditions. The propulsion system was represented in sufficient detail so that 
interactions between aircraft motions and the propulsion system dynamics could 
be investigated. 

Six degree -of- freedom aircraft motions together with the three-axis sta- 
bility augmentation system were represented. The mixed compression inlets and 
their controls were represented in the started mode for a range of flow con- 
ditions up to the inlet unstart boundary. Effects of inlet moving geometry 
on aircraft forces and movements as well as effects of aircraft motions on the 
inlet behavior were simulated. The engines, which are straight turbojets, 
were represented in the afterburning mode, with effects of changes in aircraft 
flight conditions included. The simulation was capable of operating in real 
time. 


INTRODUCTION 


One portion of the overall NASA YF-12 research program is directed toward 
investigations of interactions between aircraft motions and elements of the 
propulsion system during high-altitude supersonic cruise conditions. The 
program includes flight test, wind tunnel, and control investigations. These 
interactions become of increased importance for aircraft with mixed- compression 
inlets such as those used on the YF-12. While the mixed-compression inlet is 
the most efficient type for Mach numbers greater than about 2.2, the inlet is 
quite sensitive to aircraft motions and the resulting flow changes. Further- 
more, precise control of its moving geometry is required to achieve this design 
efficiency. This report describes a digital computer simulation which was 
prepared for the investigation of these interaction problems. The simulation 
was intended to be used for two purposes. The first was for comparisons of 
predicted with flight test results, while the second was for developing con- 
trol systems which would reduce the interactions to suitably small levels. 
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The simulation was developed for overall representations of the aircraft 
motions and propulsion system dynamics for cruise conditions at high altitude 
(above 15.24 km) and supersonic Mach numbers (greater than about 2.4). Six 
degree- of-freedom aircraft motions together with the three-axis stability aug- 
mentation system were represented. Each inlet, operating in the started mode, 
was represented in sufficient detail so that effects of aircraft motions on 
controlled inlet performance and effects of moving inlet geometry on aircraft 
motions were included. Each engine was represented in the afterburning mode 
in a linearized form together with interactions with inlet and aircraft flight 
condition changes. 

The parameter variations for the aircraft motions and engine representa- 
tions over the flight range of interest were sufficiently moderate so that 
they could have been stored in the computer over the entire flight range of 
interest. However, the generation of pertinent inlet flow parameters presented 
a special problem since the representation of multivariable nonlinear functions 
was needed. 

The use of a tabular form for these functions over the complete range of 
interest would have required an excessive amount of storage and computation 
time. In order to facilitate more efficient computer utilization, the simu- 
lation was designed to operate only in restricted ranges near selected nominal 
conditions. Nonlinear functions, principally ones describing inlet parameters, 
were represented as multivariable power series expansions. These expansions, 
then depended upon the nominal condition selected. A separate procedure was 
used to generate the expansion coefficients for selected nominal conditions 
from available wind tunnel and/or flight data. 


NOTATION 


A. 

3 

a 

b 

c 

C D 

C L 

Cl 

C 

m 

C 

n 


incremental engine primary exhaust nozzle area 
speed of sound 
wing span 


reference wing mean aerodynamic chord 

P 

V 

L 


drag coefficient, 
lift coefficient. 


q S 
n v 


^ rolling moment 
rolling-moment coefficient, 

. ^ . . pitching moment 

pitching-moment coefficient, -“g 

q v 

. . yawing moment 

yawing-moment coefficient, 
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c 

y 

. , c ££■■*. side force 

side-force coefficient, — = 

q s 

M v 

D 

aerodynamic drag 

V > 

nonlinear function contributing to inlet pressure recovery 
calculation 

V ) 

nonlinear function contributing to bypass signal pressure 
calculation 

f c c ) 

a 

nonlinear function contributing to airflow unstart boundary 
calculation 

f c < i 
s 

nonlinear function contributing to spike unstart boundary 
calculation 

g 

acceleration due to gravity 

h 

altitude 

1 , 1,1 
xx’ yy* zz 

rolling, pitching and yawing moments of inertia about body 
reference axes 

I 

xz 

product of inertia about body reference axes 

L 

aerodynamic lift 

M 

Mach number 

m 

aircraft mass 

N 

incremental engine rotor speed 

n 

z 

normal acceleration at the eg 

p»q> r 

rolling, pitching and yawing angular velocities about 
body reference axes 

p s /p t 

m 

ratio of bypass static pressure signal to free-stream total 
pressure 

p t 2 /p t 0 

ratio of total pressure at inlet diffuser exit to free-stream 
total pressure 

p t/ p t 0 

ratio of incremental burner total pressure to free-stream 
total pressure 

% 

dynamic pressure 

s 

wing area 
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W fab 

W fpb 

X,Y,Z 


x 

s 
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<S 

a 


Aw . 

1 

c 


Ax 



6 

e 


6 

r 


6 

sp 


0 


Laplace transform variable 
time constant 
total temperature 
nominal forward velocity- 

incremental airflow demanded by the engine corrected to 
diffuser exit conditions 

" e ^ T t 2 ' /T std 

PtAtd 

incremental afterburner fuel flow 
incremental primary burner fuel flow 
forces along the x,y,z body axes 

incremental shock wave position from unstart location 
incremental angle of attack from a 

w o 

angle of attack relative to wing reference plane 
sideslip angle 

total antisymmetrical elevon deflection 

inlet airflow connected to diffuser exit conditions, 

W i 

denotes incremental value from the nominal condition, x-xo, 
for the arbitrary quantity x unless otherwise indicated. 

position of actuator controlling bypass exit flow 

symmetrical elevon deflection 

vertical tail deflection 

inlet spike position 

incremental engine thrust 

pitch angle 
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p 


mass density of air 
roll angle 


♦ 


Subscripts 

a 

e 

m 

o 

s 

sas 

std 

2 

4 


total antisymmetrical values of left and right side 
propulsion system parameters 

external input 

measured value 

nominal value 

total incremental symmetrical value of left and right side 
propulsion system parameters 

output from stability augmentation system 

standard value at sea level 

value at compressor inlet 

value at primary burner exit 


Superscripts 

£ 

r 

s 


left side 
right side 

stability axis system 


Aerodynamic Derivatives 

Generally, derivatives of coefficients are with respect to the indicated 


dC 


subscript; e.g. C 


n 


n„ 


d|3 ‘ 


An exception to this notation is used for the 


aerodynamic damping derivatives , , Cj , . These terms are deriv- 


dC, 


atives with respect to the subscripts times e.g. C 


'1 


d(rb/2V) * 
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DESCRIPTION OF AIRCRAFT 


The YF-12 is a delta-wing twin-engine aircraft which is capable of 
cruising at Mach numbers greater than 3 and at altitudes greater than 
24.38 km. A three-view drawing of the airplane is shown in figure 1. 

The aerodynamic controls are as follows. Two elevons on each wing, one 
inboard and one outboard of each nacelle, operate together symmetrically and 
antisymmetrically with the pair on the opposite wing to provide longitudinal 
and lateral control. Two nacelle mounted all-movable vertical tails provide 
directional control. 

The aircraft has two axisymmetric variable geometry, mixed compression 
inlets, Each inlet has a translating centerbody (spike), which varies the 
inlet contours with flight condition and a forward bypass door which diverts 
flow from the engine and dumps it overboard. The forward bypass door is 
continuously modulated through use of a closed loop control system to main- 
tain the inlet airflow at a desired condition and to prevent inlet unstarts. 
Mach number as well as angle of attack and sideslip angle signals are sup- 
plied to the spike and bypass systems. Other flow is diverted from the inlet 
by a fixed bleed system on the centerbody which exhausts the flow overboard, 
a fixed bleed on the cowl, and an aft bypass door. The latter two flows supply 
secondary air for engine cooling. The aft bypass doors are manually scheduled 
by preselected increments with flight condition. In the event of an inlet 
unstart, control signals are scheduled to adjust the inlet spike and forward 
bypass doors to restart the inlet flow and then return the inlet to its 
previous position. 

The engine is a Pratt and Whitney YJ58 single-spool axial flow turbojet 
with afterburner. Provision is made to bypass a portion of the air from the 
compressor. The bypassed air remixes with the primary flow at the entrance 
to the afterburner. Engine rotor speed is controlled by closed loop adjust- 
ment of the exhaust nozzle which varies back pressure on the turbine. The 
flow exits through a converging-diverging ejector nozzle with adjustable 
trailing edge flaps. Afterburner fuel flow is an open loop function of power 
lever angle, burner pressure, and compressor inlet temperature. Primary 
burner fuel flow depends upon the same quantities plus rotor speed. 


DESCRIPTION OF SIMULATION 


The purpose of the simulation was to represent both airframe motion and 
propulsion system dynamics in sufficient detail so that effects of parameter 
variations and interactions between the systems could be investigated. The 
simulation was intended for use both for identification with flight test 
results and for control system studies. Since a simulation which would 
operate at least as fast as real time was desired, a careful balance was 
necessary between depth of representation of the various components and com- 
putational speed. The simulation was intended to operate in the aircraft 
motion frequency range, which included the phugoid mode, from about 0.005 to 
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to 5 hertz. Hence higher frequency effects associated with both aircraft 
motions and the propulsion system were neglected. Simplifications in the 
representation of several inlet nonlinear functions, obtained from wind tunnel 
data, were also necessary. Rather than attempt to represent these functions 
of several variables in a tabulated form for the complete range of interest, 
power series expansions, valid for a restricted range, were made relative to 
selected nominal flight conditions. A procedure for obtaining these functions 
is discussed in a subsequent section and Appendix A. It was intended that the 
expansion coefficients be obtained through use of a separate computer program 
which would curve fit the original data and obtain the coefficients in a form 
for direct use in the simulation. 

The equations used for the simulation are described in detail in the 
following sections of this report. Descriptions are given of the aircraft, 
inlet, and engine representations, as well as their respective control systems. 
Then the numerical integration procedures used are presented. Finally, a 
listing of the computer program is provided in Appendix B. A sketch of the 
subsystems represented together with the interconnecting variables for these 
systems are shown in figure 2. The numerical data used for the simulation 
are not included in this report. 


Airframe Representation 

The equations used to obtain six degree-of-freedom rigid-body aircraft 
motions are given in this section. Effects of inlet and engine forces and 
moments are included. The representation for the aircraft three-axis stability 
augmentation system and the control surface dynamics are also described. An 
outline of the computations performed is shown in figure 3. 

Aircraft equations of motion — The following assumptions were made in 
determining the airframe equations of motion. 

1. The airframe is assumed to be a rigid body. 

2. The mass of the airplane is assumed to be constant for the portion 
the flight to be analyzed. 

3. The earth is assumed to be flat, nonrotating, and fixed in inertial 
space. 

4. The aircraft has a vertical plane of symmetry. 

5. Small angle approximations apply. 

Two moving- axis systems were used with origins at the center of gravity 
of the aircraft. The first is a stability axis system. In the definition 
used here, the x s and z s axes (positive forward and downward respectively) 
are located in the aircraft plane of symmetry and are aligned parallel and 
perpendicular to the instantaneous relative wind projected into the plane of 
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symmetry. The y s axis is perpendicular to the plane of symmetry (positive to 
the right). The second axis system is a body axis system which is fixed to 
the aircraft. The x and y axes are fixed in a wing reference plane (denoted w) 
which is aligned with the principal wing surface, whereas the z axis is per- 
pendicular to it. The transformation from one axis system to the other is 
made by a rotation of the angle, a w , about the common y axis. Aerodynamic 
stability derivatives and force and moment coefficients were available in the 
stability axis system. These quantities were transformed to the body axis 
system and then integrated with the result that motions as measured in flight 
would be directly represented. Moreover, body motions to be used as SAS sig- 
nals would also be directly obtained. Since effects of changes in M on sta- 
bility derivatives are small for the range of interest, they have been 
neglected. 

As part of the initialization procedure, lateral-directional stability 
derivatives were converted from the stability axis direction to the body axis 
direction (a rotation of a^) . The effect of any subsequent changes in a on 

the lateral- directional stability derivatives was neglected. For these static 
derivatives, the conversion from stability to body axes was as follows: 

C = C S cos a + sin a (1) 

n. n. w Z. w 





sin a 

w 


o 


( 2 ) 


where the subscript, i, denotes a static derivative such as 3, 6 a , etc. For 
the damping derivatives, the equations for resolving both. the moment coeffi- 
cients and the angular rates from stability to body axes are: 




„s . 2 
C s a 
n w 

P o 


a 

w 

o 




sin a 

w 


o 


cos 




+ 



(4) 


(5) 
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( 6 ) 


C. = C S sin 2 a +■ cos 2 a 


'SL n 
P r 


w i w 

op o 

' \ 

C + C„ sin a cos a 
n £ / w w 

, p r/ o o 


where the subscript, s, denotes stability axes and no superscript denotes 
body axes. The body axis derivatives were used in all subsequent calculations. 


Drag and lift forces were initially calculated in the stability axis 
system. 


— = C®(a) + cf. A6 + <5, + C® <5 

m m D’ k ‘ ; D. e D- bp D. sp 

6 6, *s <5, *s 

L e bt> bp J 


D. e bp D. sp 

6 6, *s <5, *s 

e bp bp 


— = ilL. [c? ( a ) + cf AS + C? 5. + C? 6 

m m L l L„ e L. bp L. sp 

6 <S. *s o r s 

L e bp sp 


C -S 

* 2V C L <> 

q J 


These forces were subsequently converted to body axis directions in the equa- 
tions of motion. 


The following equations describe the aircraft linear and angular accel- 
erations relative to the body axes. 


i» + \ v q ■ -g A ® + * 6 6 t 


D S L S 

— + — a + X 
m m w o 


q v S 

v + Vr - g sin (j) - a Vp = 


■ C 3 + C 5 i + v, 

m \ ^ r 2Vm 


% Sb 

o 


w - Vq = Z x 6 - a - — + Z 

“t *s m w m c 


I q Sb / ' 

p _ _i- r = j (C 0 8 + C. a8 + C 6 + C„ 6 + C„ 6. + C„ 6 

*xx *xx 1*8 *aB *6 3 *6 r \ b P a *6 S P a 

\ a r bp sp / 

q Sb 2 

Mr i V 


C i * C «. P * L S 
r p / t a 
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q v Sc 


q Sc n 

V(c.i*c 1 ,)-- f 2-c 1 (0) 


q 2 VI l“m. 

yy \ a 


yy 


XZ 


q Sc 
n v 

q Sb 


C (a) + C AS + C 

V J »T| ^ 1 


m ' ' m t e bp g 


bp 




P * 


zz 


zz 


C$ + C 6 + C S + C 6, +C <5 

n e n 6 3 n 6 r V bp a n s sp aj 
p a r bp sp 

Sb 2 


2VI 


zz 


: + c pi + n, 6 

n n r / cit 


t a 


( 13 ) 


(14) 


Total variables were defined from incremental and reference condition 
variables as follows: 


Propulsion System Variables: 



6 










r 


where i denotes a spike, bypass, or thrust variable. 


( 15 ) 


Longitudinal Variables: 


0=0 + A0 

o 


6 = 5 + A<5 I 

e e o e 1 


(16) 


The horizontal and vertical forces for the reference conditions in equa- 
tions (9) and (11) were calculated from the drag and lift terms by the fol- 
lowing equations. 
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( 17 ) 


% S 
x = -~sr 

o m 
q v S 

z o - ~m (V 0) * C D (0) \) I 

Note that the transformation from stability to body axes for the lift and 
drag longitudinal equations was approximated differently than for the yawing 
and rolling moment lateral-directional equations. Since the effect of angle 
of attack changes were generally greater for the longitudinal equations, the 
lift and drag equations were continuously transformed from stability axes to 
body axes (eq. (9) and (11)) using the instantaneous values of angle of attack 
rather than only the nominal value. However, small angle approximations were 
used. 


(' 


c D (°) 


C T (0)a 
L w 


As indicated in previous equations, only aerodynamic nonlinearities due 
to angle of attack were included. These quantities were expressed in power 
series form as follows: 


= d + C® a + d ct 2 
D- D w D 2 w 
fa or 


dja) = d + d L a w 

fa 


C(a)-C +C a +C a 
m y nu m w mo w 
fa a 


(18) 

(19) 

( 20 ) 


where the subscript f refers to the value of the quantity at a w = 0. 
Note that from the definitions used 


a 


w 


a + a 
w 

o 


Dynamic pressure was calculated from the following approximate equation 


q v = \ p(h)(V + u) 2 


( 21 ) 


For the nominal flight condition, the equation becomes 

1v = 1 “ V 2 

o 
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The mass density, p(h), was obtained from a tabulation of the 1962 atmosphere 
(ref. 1). The variable dynamic pressure, q y , was replaced by the nominal 
value, q v , for the smaller quantities such as damping terms in equations (9) 
to (14). 0 

Incremental Mach number was calculated from the following approximate 
equation. 

AM = u/a Q (22) 

Small angle assumptions were used for the Euler angle transformations, 
flow angles, and vertical velocity equations. 


4* ■ P 

(23) 

e = q 

(24) 

w 

a - v 

(25) 

ft _ V 

e - 7 

(26) 

h = v (e - - v<t> 

(27) 


Definition of -propulsion system coefficients— Since the simulation was 
concerned with interactions between the propulsion system and aircraft motions, 
effects of both inlet and engine exit flows on the airframe forces and moments 
were included. In equations (9)- (14), these effects appear as the bypass, 6bp* 
and spike position, 6 sp , teims from the inlet, and the thrust, <5 t , from the 
engine. The following interpretation of these propulsion system terms was 
used. The external forces and moments on the airframe due to an element of 
the propulsion system can be developed from external flow changes associated 
with that element. The resultant forces correspond to changes in momentum and 
energy of the discharged flow from that present in the undisturbed flow 
upstream of the aircraft. Changes in the internal flow are represented in the 
propulsion portion of the simulation. For example, a change in bypass position 
has two effects. The first is a change in external forces which result from 
energy losses and direction changes in the flow discharged from the bypass exit 
This effect appears as the bypass terms in equations (9)- (14). The second is 
a change in the internal inlet flow, which is calculated in the propulsion 
part of the simulation. This change is transmitted to the engine, and results 
in a second external force change due to engine thrust. Similarly, a change 
in spike position produces a change in the flow at the entrance to the inlet 
and a resulting change in external forces, as well as a change in the internal 
inlet flow to the engine which results in a thrust change. The external 
effects of the spike are relatively small in comparison with the bypass effects 
While the external forces due to the bypass and spike position vary somewhat 
with angle of attack and Mach number, a linearized representation suitable 
for a restricted range of these variables was used with the simulation. 
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Note that when effects of inlet geometry changes on overall forces are 
obtained from wind tunnel results, care is needed to insure that the measure- 
ments are consistent with the definitions used for the simulation. That is, 
internal inlet forces need to be measured and subtracted from the overall wind 
tunnel values so that only external forces remain. Two additional factors 
must be considered in the application of wind tunnel results to the represen- 
tation of flight results. The first is that the internal pressure recovery, 
and hence the energy of the discharged inlet flow will generally be less for 
the model than for the aircraft. The second is that the angle of the dis- 
charged inlet flow for the wind tunnel model may be different from that on the 
aircraft. Hence approximate corrections may be needed for the wind tunnel 
results to account for these differences in the aircraft representation. 


SAS and control surface dynamics— The pitch, roll, and yaw SAS dynamics 
for the aircraft were represented for the high-altitude supersonic cruise 
conditions of interest. The system provides basic damping of the aircraft 
motions and utilizes primarily angular rate sensing. In addition the yaw SAS 
is designed to provide not only airframe damping, but also stabilization about 
the vertical axis during an inlet unstart. The yaw SAS uses a lateral accel- 
eration sensor at the nose as well as a yaw rate sensor at the c.g. Effects 

of sensor dynamics and system gain changing with flight condition were 

neglected in the simulation. The equations used are as follows: 


p 

*sas 


r 

sas 


K p 
PP 

1 + v 

T ri K rr rs + 
1 + T ri s 


K n 
m y 
7 n 

1 + T s 
r2 


(28) 

(29) 


where n is the lateral acceleration at the nose 


n 


q = q 
sas n 


rs£ : v s) 

|_ <> * 


yi : 

(1 * T q3 s)U ♦ T^s) 


(30) 


Limits were included on the total p S as» 'Isas* an d r sas f erm s, and separately 
on the second order term in the q sas equation. 

The servo actuator dynamics for the aerodynamic control surfaces were 
approximated by second order equations. Provision was made in the simulation 
for SAS inputs, initial conditions, and external inputs. Effects of aero- 
dynamic forces on the surfaces were neglected. The equations were expressed 
in the following form. 
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Aileron: 


(31) 


Rudder: 


Elevon : 



2? a . 

6 +6 = p + p 

a) a a r sas r e 

n 

a 

25 r • 

6 + 8 = r + r 

cj r r sas e 

n r 


2r, 

— — 6 + A6 - q + q 

<ii e e ^sas n e 

n e 


(32) 


(33) 


Inlet Representation 

This section describes the simulation of the inlet together with the 
interconnections with other portions of the system. The first part provides 
an overall description of the inlet and the variables used to represent it. 

The succeeding parts describe several portions of the simulation in more 
detail. 

Description of inlet system- The inlet for the YF-12 aircraft is of the 
mixed compression type with a translating centerbody used for geometry changes. 
Although the inlet is basically axisymmetric , body interference effects, 
together with some geometry changes to adjust for these effects, result in 
characteristics which are somewhat unsymmetric in the presence of sideslip 
and angle of attack variations. Forward and aft bypass doors are provided for 
control of the inlet airflow. The forward door continuously modulates the 
flow whereas the aft door is adjusted manually to one of three positions as a 
function of Mach number and serves only as a coarse flow adjustment. Effects 
of only the forward door were included in the simulation. The inlet also has 
a system of bleeds on both the cowl and centerbody. Airflow from these bleeds 
is either dumped overboard or is ducted to the engine to provide cooling. 
Provision for variations in these secondary flows was not included in the 
simulation. Effects of inlet flow distortion and turbulence, which are of 
importance for inlet-engine interactions, were also neglected. 

A block diagram of the inlet simulation is shown in figure 4. The repre- 
sentation for each inlet is the same except that the asymmetrical 0 functions 
for one side are the mirror images of the 6 functions for the other side. The 
simulation is intended to represent conditions for the started inlet up to the 
unstart boundary. Inlet flow dynamics for the diffuser were simplified to a 
first order lag term since only motions in the aircraft frequency range were 
of interest. Hence the principal dynamics represented were those due to con- 
trol system components. The basic parameter used to represent the state of 
the inlet diffuser was inlet corrected airflow, Wj^; i.e. diffuser airflow 

corrected to diffuser exit conditions. This parameter is a direct function of 
inlet exit geometry for the frequency range of interest. The exit changes 
result from variations in the forward bypass opening and from changes in the 
airflow demanded by the engine. 
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Three other inlet parameters were determined - (1) pressure recovery, 

(2) bypass signal pressure ratio, and (3) terminal shockwave position. Inlet 
pressure recovery is a parameter affecting engine performance and is a primary 
input to the engine. The bypass signal pressure ratio is needed for the 
representation of the bypass control system and terminal shock-wave position 
represents the margin of the inlet from the unstart condition. A major problem 
in the simulation was the incorporation of wind tunnel results to provide an 
adequate representation of these three quantities. Wind tunnel data included 
in a Lockheed Aircraft Corp. report of a detailed digital computer simulation 
of the inlet indicated that the quantities are nonlinear functions of five 
variables. The two variables which have primary effects are the inlet vari- 
ables: corrected inlet airflow and spike position. Three more variables, 

which have a secondary effect, are those describing entrance conditions to 
the inlet: angle of attack, sideslip angle, and Mach number. A complete 

tabular representation of functions of this many variables, even over 
restricted ranges, would result in a rather cumbersome and slow computational 
procedure. A simpler approach for the purpose of generating time histories 
was to first curve-fit the multi-dimensional wind tunnel data by a power 
series expansion in several variables. Each function was then determined from 
the series as needed during the computation of the time history. 

Inlet pressure representcctvon— The coefficients of the series repre- 
senting the pressure recovery and bypass signal pressure functions were 
obtained through use of a least square criterion. Care was necessary in 
selecting the range to be curve fitted and the resulting number of coefficients 
to be used in order to keep the representations as simple as possible. As 
previously mentioned, these pressures were represented as functions of two 
primary variables, the inlet variables, w i c > and ^sp» an< ^ three secondary 

variables, a, g, and M. The functional representation was expressed as fol- 
lows. For a nominal condition of the three secondary variables, an expansion 
was made of the two primary inlet variables, Wj c , and ^ S p* Expansions of the 

primary variables were also made with each of the secondary variables. Note 
that this procedure neglects terms containing crossproducts between the sec- 
ondary variables. In these expressions, the expansions for the inlet vari- 
ables, Wi and 6 , were made relative to convenient average inlet conditions, 

C 

while the expansions for the secondary variables were made relative to the 
particular nominal flight condition. The former choice was made in order to 
facilitate adjustment of the inlet operating condition. Hence, for a set of 
nominal values of u and M (g = 0) , the pressure recovery and bypass signal 
pressures were determined by the following expressions (fig. 4) : 

— = f /w. , 6 'j + 

P t a V V S *7 
o 

+ f | w. , 6 , 

a V i c ’ sp’ 


f a (“I/ V “) 

b) . f a . m) 


(34) 
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( 35 ) 



The series expansions 



for each of the f functions were expressed as 


) iricLX in a. a. iii ua 

'EEE 


7 max max 


c iik"i c % 


j x k 


i=0 j=0 k=l 


(36) 


where the quantity, X, represents a, 3, or M. The upper limits used, which 
determine the number of coefficients required for each pressure function, are 
listed in Table 1. In the selection of the coefficients, Cjjk, an attempt 
was made to provide a balance between accuracy required and the need to main- 
tain computational time as short as possible. Interpolation of the functions 
was also necessary in order to obtain expansions about the desired nominal 
flight conditions. The procedures used for determining these coefficients 
are described in Appendix A. 

As previously stated the input variables were expanded relative to con- 
venient average inlet conditions. Hence, for a particular reference condition, 
steady state corrections were necessary to set the nonlinear functions at the 
correct values and to interface with linear portions of the simulation. The 
corrections provided a means for adjusting the inlet to different reference 
conditions. The quantities, 6 S p o and wi CQ (fig. 4), were used to define the 
inlet nominal operating conditions and, hence, to calculate nominal values of 
the nonlinear pressure and unstart functions. The quantity, (Ps/Pt m ) Q » was 

calculated to null the signal error for the inlet nominal conditions so that 
a proper interface with the incremental bypass position, A^bp* would be 

obtained. The quantity, ^p t2 /Pt 0 ) 0 was ca l cu ^ a ^ e< ^ to nu ^ th e pressure 

recovery signal to the linearized engine for nominal conditions. Note that 
these corrections are distinct from initial condition selections. 


Shockwave position and unstavt boundaries— kn additional calculation was 
made to obtain the inlet unstart boundaries and shockwave position as a func- 
tion of distance from the unstart condition. The calculation was based on a 
power series expansion of wind tunnel results (contained in a Lockheed Air- 
craft Corp. report) with respect to a nominal flight condition. The form for 
the inlet unstart boundary to be represented as a function of inlet parameters 
(with fixed a, 3 and M) is shown in sketch a. The lower boundary in the sketch 
represents the effect of restricting the inlet exit flow so that unstart occurs 
and will be called the airflow boundary. The left hand boundary represents 
the effect of over contracting the area ratio at the inlet entrance with the 
spike so that unstart occurs, and will be called the spike boundary. These 


16 


A- 4840 




boundaries are affected primarily by changes in M and also by variations in 
a- and $. The data for the unstaxt boundaries were obtained from steady-state 
wind tunnel tests. Because of the relatively low frequency of the disturbing 
quantities relative to the basic inlet dynamics, dynamic effects on unstart 
were not considered significant. 


An examination of available wind tunnel data resulted in the following 
functional forms for the unstart boundaries. The relation of each function 
to the rest of the inlet simulation is indicated in figure 4. In curve fit- 
ting the data, effects of a, 8, and M were represented separately since 
effects of simultaneous variations in these quantities were felt to be of 
higher order. The basic airflow and spike boundaries, with M corrections 
added, were expressed as follows: 

Airflow Boundary, f £ (6 \ f Q U , mV 


For 6 < 

s Pl 


6 <6 
sp d sp 2 


For 


5 > 

sp d 


sp 2 



where 6 = <$ + K. AM 

sp d sp h 2 


f M 1 = -K. AM + K. 6 AM 

C a\ Sp / "15 he sp 


Spike boundary. 





6 


sp 3 


(37) 


(38) 


(39) 
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All of the quantities in equations (37) and (39) with numbered subscripts 
denote constants. The M corrections provide both slope and displacement 
changes to the airflow unstart boundary. 

Separate a and 0 functions were applied to obtain displacement corrections 
for the airflow and spike unstart boundaries. These functions were of the 
following form: 


Angle of attack correction: 


Sideslip angle corrections: 
For 0 < 0, 

i 

c 

For 0 > 0, 


f c («) - « i26 a » W 

a 


f («) - K. a + K. a' 

C 128 129 

S 


f- CM - K i2 ,6 2 ; f c cm = x 124 e* 


f. cs) - f. CM - k. 25 b 


122 


i 


(40) 


(41) 


As indicated in figure 4, the unstart airflow boundary, w i Cu > is the sum 
of the airflow boundary nonlinear functions. 


w. = f (t \ + f I*. A * f (a) * f (0) 
1 c u C a\ s Pd/ C a\ S P ) C a 


(42) 


Similarly, the unstart spike boundary is the sum of the spike boundary 
functions. 


6 


sp 


u 



+ f Ca) + £ (0) 

c s c s 


(43) 


The shockwave position was curve-fitted from wind tunnel data as a non 
linear function, f d (w i( J , of incremental airflow from the unstart boundary. 


where 



(44) 
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The function was assumed to be independent of inlet and flight conditions for 
the operating range of interest. An inlet unstart, which could result from 
either the airflow or spike boundaries being exceeded, was indicated by either 
of the following equations: 


* s 6 0 , 


or 




(45) 


The simulation was terminated when an unstart occurred since unstart conditions 
were not represented. 

Inlet aontrol system— The inlet control system positions the spike as a 
function of flight condition, and bypass flow as a function of flight condition 
and engine airflow demand. The objective is to maintain the inlet airflow at 
as high a performance condition, i.e. high pressure recovery and minimum flow 
distortion, as possible while still preventing an inlet unstart. The spike 
position is determined by an open loop command function which is computed from 
flight condition measurements. Bypass flow is adjusted by closed loop control 
with the forward bypass doors. The bypass signal is a pressure signal indica- 
tive of the position of the terminal shockwave in the inlet. It is formed from 
a combination of several static pressure sources which are summed and divided 
by a signal which is an approximate measurement of free stream total pressure. 
The use of this latter measurement reduces effects of dynamic pressure changes 
on the control signal. The signal is represented by a single variable, 
p s /pt m > in the simulation. The desired pressure signal is a function of 

flight condition measurements. 

The inlet control command signals and bypass loop are included in fig- 
ure 4. The flight condition measurements, a, B, M, and n z were used as open 
loop commands to both the spike and bypass controls. Dynamics associated with 
these measurements were represented by the time constants, T a , Tg, and T^. 

The effect of a, 3, M, and n z command signals to the spike were represented 
by the functions f S p(“ m )> K il2 , f(B), K il9 , and f sp (n z ) respectively. Simi- 
larly command signal corrections used for the bypass pressure ratio setpoint 
were f b p(otM), K ilQ , f (B) , and f^p(n z ). The command functions were expressed 

as power series expansions as follows: 


Spike commands: 


where 


4 £ s p (\) ' fs e(°”J ' ^(“0 
/ \ 3 

f ( < a | = \ a . | a - 5 . 0 1 1 

S P\V 4-f 301 W m 

1=2 


(46) 
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f (m J = K. .AM 
sp\ m/ il9 m 
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For |n z | > £i> 
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For |n z | < , 
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f (n 1=0 
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where is a constant. 
Bypass comman ds : 
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m 
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m 
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( 47 ) 


(48) 
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Spike and bypass command: 




For 6 m < 0, 

f (0 

For 6 m > °» 

f (0 


£ a 6i 8 m 


1 = 2 




i = 2 


1 

SiTn 


( 51 ) 


Note that the command functions are expressed as incremental values from nom- 
inal conditions. 


The bypass control loop consisted of the bypass control dynamics and the 
previously described diffuser dynamics and signal pressure function (fig. 4). 
The bypass actuator dynamics and filtering were simplified to the second order 
form shown. Effects of M changes on the bypass control gain were neglected. 


Engine Representation 

The simulation represented the dynamic response of the engine operating 
in the afterburning mode and a block diagram is shown in figure 5. The simu- 
lation for each engine was identical. The response was linearized relative 
to a selected nominal condition. Hence, all engine states shown in figure 5 
are incremental values from the nominal condition and the A symbol denoting 
increments has been omitted. The primary purpose of the engine representation 
was to provide input- output relationships for those variables which directly 
interacted with other portions of the aircraft systems. However, sufficient 
depth of the representation was considered necessary to also provide a basis 
for comparisons with subsequent flight-test results. The principal relation- 
ships represented consisted of the dynamic response of engine thrust, com- 
pressor airflow, and fuel flows to throttle position and inlet pressure 
recovery. Response to the latteT input was included since it represented the 
principal interconnecting parameter from the inlet to the engine. The com- 
pressor corrected airflow, w e , response was needed since this quantity was the 
the principal interconnecting c parameter from the engine to the inlet. 

Effects of variations in flight condition, i.e. Mach number, M, and 
dynamic pressure, q v , on the engine response were also represented. These 
quantities were input to the engine simulation at the pressure disturbance 
location with the gain functions, Ki 23 and fj(M) based on calculated thrust 
changes. The gains were determined from steady-state variations in calculated 
thrust due to each input parameter with the other parameter held constant. 

The quantity, q v , produces an engine dynamic response which is similar to a 
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p t ^/p t input since a compressor inlet pressure change occurs in both cases. 

Hence, these two inputs were added at the same point for the simulation. How- 
ever, the M disturbance produces engine dynamic responses somewhat different 
from’pt /p t since it introduces compressor inlet temperature changes. Since 

the frequency content of the M disturbance is relatively low in comparison 
with the engine dynamics, the representation of the difference in the dynamic 
response due to this disturbance from a Pt 2 /Pt 0 disturbance was not considered 

necessary. Hence, the M disturbance was introduced in the same manner as the 
pressure disturbances. As previously mentioned, the gains from the q v and 
M inputs were determined from the variations in steady-state thrust only. 

Hence, although qualitatively the same, the resulting changes in fuel-flow and 
burner pressure due to M inputs were somewhat in error. Note, however, that 
these engine variables are of less importance than others since they are not 
fed back to other portions of the simulation. 

Corrections due to M inputs were also applied to the engine compressor 
airflow variable by means of the gain Ki 20 . The gain was estimated through 
use of a linearization of steady-state compressor operating curves and the 
corrected airflow relationship. An additional consideration was that a change 
in compressor inlet temperature due to M would result in a slight change in 
trimmed rotor speed. This effect was neglected in the rotor speed calculation 
but was included in the selection of the gain for the engine airflow calculation. 

Once the form of the engine simulation was established, numerical values 
of the dynamic coefficients were obtained by matching time histories computed 
from this simulation with results computed from a more extensive simulation of 
the engine which included nonlinear representations of all major engine com- 
ponents. The engine simulation is described in a report by the Pratt § Whitney 
Div. of United Aircraft Corp. For a desired nominal operating condition, 
responses of installed thrust, fuel flow, burner pressure, and rotor speed to 
throttle and pressure recovery inputs were used for the comparison. 


Numerical Integration Techniques 

The methods used for numerical integration of the dynamic equations 
previously described are given in this section. In selecting the integration 
method to be used for achieving a certain accuracy level, careful attention 
must be given to tradeoffs in computational time, levels of complexity in the 
method itself, and the maximum allowable size of the integration step. This 
tradeoff will also depend on the degree of complexity of the system repre- 
sentation including the nonlinear relationships present. 

Two basic integration methods, which were intended for use with fixed 
step size integration, were applied to different portions of the simulation. 
The first was an equation of motion method for integration of a set of first 
order differential equations. This method was used to integrate the aircraft 
equations of motion. The second method involved the determination of the 
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response of linear portions of the system to a step input during each integra- 
tion interval and is called the z-transform method. This method was used for 
the remaining portions of the system - SAS, inlet, engine, and their controls. 

For the equation of motion method, the following second order Adams- 
Bashford predictor equation (e.g. ref. 2) was used to obtain the value of the 
vector state Y at time t^ +1 from known values at times t k and t^. 

*(V.) - Y ( l k) * [ 3Y '( t k) - Y '( t k-i)] (T/2) (52) 

where the continuous 1st order differential equation is given by the vector 
form 


Y* =* f (t,Y) 

The method involved prediction by means of a polynomial curve fit of present 
and past values of each state variable to be integrated. While a higher order 
integration algorithm involving more past terms in Y* could have been used, 
some computational experimentation indicated that an overall improvement 
would not be obtained. In order to maintain a certain accuracy level, the 
computation interval could not be increased sufficiently to compensate for the 
additional time required to compute the higher order form. A basic problem 
in using this integration method occurs when a large dynamic range of equations 
is present. If the step size selected is too large for the higher frequency 
dynamics, the calculation can go unstable since the derivative for the high 
frequency equation is assumed constant for too large an increment of time. 

One method for improving computational accuracy for the higher frequency 
portions of the system is to integrate those equations in closed form with the 
assumption of a known input over the integration interval. The z-transform 
method can be used to find the output of a linear dynamic element to a selected 
input over a sampling interval, T. This approach is particularly useful for 
a system represented in block diagram form. An example of the application of 
the z-transform method to numerical integration is included in reference 3. 

Note that the method can still be used for a nonlinear system as long as the 
system can be separated into zero memory nonlinearities and linear dynamics. 

The simplest form of input to use is a step that is assumed constant over the 
integration interval. However, the input actually changes during the interval. 
If the value of the input at the beginning of the interval is used for the 
entire interval, a time delay is introduced which can result in an unstable 
calculation for a system with feedback loops. One way to include the effect 
of the changing input, with values known at the beginning and end of the 
interval, would be to determine the response of the linear element to a ramp 
input. A simpler approach, which was actually used, is to employ the equations 
for the step response but to average the input over the sampling interval. 

This averaging was done as follows. For locations where the input was 
known at both the beginning and end of the sampling interval, the average 
value of the input was directly calculated by the following equation. 
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r av = I [ r (*k) * r (Vl)_ 


(53) 


This computation was possible for instance when integrating a state with an 
input which was previously integrated through use of the equations of motion 
method. However, the value of the input at the end of the sampling interval 
may not always be available. For instance at feedback junctions, the updated 
value of the returned quantity generally is not known. A typical situation 
is shown in sketch (b). 



Sketch (b) 

The input, r, was assumed known up to and including the time, States in 

the feedback loop dynamics including the output, c, were known only up to the 
time t^. In order to average both inputs to f^fs), an estimate of Cj for 
the interval t^ to t^ +1 was needed. The following predicted average equation 

was used for this situation. 

c i av = i [’'iW - c i(vi)] (s4) 

P 

Again it was felt that the use of a higher order prediction equation would not 
allow a sufficient increase in sampling time to compensate for the resulting 
increased computational time. The predicted coefficient, c, was then com- 

av p 

bined with a time averaged r coefficient to form an averaged error input to 
fi(s). The integration algorithm was then applied to obtain crtfc+iV Next 
the outputs of the remaining quantities in the loop, f(c) and l^Cs) w e re com- 
puted for the time t^+j. While the procedure could have been iterated further 
through use of the quantity, Ci^t^+i), it was not considered necessary. Simi- 
lar procedures using input averaging and prediction, where necessary, can be 
applied to more general feedback structures, and also to systems of first 
order equations. 
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For the YF12 simulation, the airframe motion dynamics were expressed in 
first order form and integrated by means of the equations of motion method 
(eq. 52). All other components, the SAS, inlet, and engine, and their con- 
trols, were integrated by the z- transform method with the components in block 
diagram form. A subprogram was used to compute z-transform coefficients which 
represented responses of first or second order transfer functions to a step 
input. Input averaging (eq. (53) and (54)) was included for cases in which 
the frequency content of the input was sufficiently high to result in some 
variation during the integration interval. The simulation operated in real 
time with a step size of 0.02 seconds. 


CONCLUDING REMARKS 


The simulation was intended for use both for the identification of param- 
eters by comparison with flight test results, and for control system studies. 

Only preliminary flight measurements of SAS-off Dutch-roll oscillations have 
been available and these measurements allow a comparison with only a portion 
of the simulation. The flight measurements indicated a reduction in Dutch-roll 

damping due to activity of the inlet control system and results from the simu- 

lation indicated similar changes. Hence, on the basis of available flight 
data, the simulation is felt to exhibit adequate sensitivity to parameter 
variations. 
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APPENDIX A 


REGRESSION ANALYSIS FOR NONLINEAR PRESSURE FUNCTIONS 


The requirements for determining the coefficients for the polynomial 
expansion for the inlet nonlinear pressure functions are outlined in this 
section. The effort involves selecting the form of the expansion and obtaining 
a least squares fit of the data for the range of interest. To date these 
computations were performed separately although a general program could be 
prepared to compute all necessary coefficients in sequence and output them in 
a punched card form suitable for use in the 8400 digital computer program. 

The necessary mathematical steps are outlined in the following paragraphs. 

The polynomial expansions of the pressure variables were obtained from 
wind tunnel data which are functions of five independent variables. Two vari- 
ables were considered to be primary ones, w^ c , 6 S p, to be fitted with greater 

accuracy, and the remaining variables, o^, 8, M, were considered secondary and 
were fitted less accurately. Polynomials of prescribed form were fitted to 
the data using a least square criterion. A factor to be accounted for in the 
curve-fitting process was that the range for which the inlet variables (wi c > 

6 sp ) were needed varied with flight condition. A double interpolation of a 
portion of the polynomial coefficients was also required since, in general, the 
expansions were made relative to nominal conditions which were different from 
the original data conditions. 


The complete polynomial expansion for either the pressure recovery or 
signal pressure function, P, was expressed in the following form by the 
substitution of equation (36) into equation (34) or (35) 



■Snax J max 


y y C..W 1 

Z—J Z—t 11 1 sp 
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Note that the expansion is made relative to nominal flight condition of 

cty, = cty 8 = 0, and M = M 0 . Also, the incremental angle of attack is defined 

as a = cty, - a WQ . As previously discussed, provision was made to expand the 

inlet variables, w ic and 6 sp , about any convenient average inlet condition. 

The upper summation limits for each pressure function are given in Table 1. 

The procedure for obtaining the coefficients was separated into the following 
steps. 

The first step was to perform a linear regression for the two primary 
variables, and 6 sp , f° r va l ues o f a w , M for which wind tunnel data were 

available. Hence, the form of the first summation in equation (Al) , denoted 
by Pi, was used to determine the regression equation for the unknown coeffi- 
cients Cij(a w , 8, M). 

p. ( v e. M - J ) - 

i=0 

where 

J = 1 • * * N^a w , 8, MJ 

The value, N^cty,, 8, m), is the number of data points used to determine 
p l(«w> M ) as a function of vx c and 6 sp for each constant (a w , 8, m) 
condition. 

The second step was to interpolate the coefficients in the a and M 
directions to selected nominal values of and M Q . Both linear (using two 

values per direction) and quadratic (using three values per direction) inter 
polations were used. The choice depended upon the variations in the coeffi- 
cients between data points in the and M directions. An available single 
direction interpolation program was successively applied as follows. 



where the subscript o refers to the selected nominal value of the M or a w 
variable. 


L c u(v e - m K (A2) 
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The third step was to perform a linear regression for the primary vari- 
ables with each of the secondary variables, a, 3, M, as indicated by the 
second, third, and fourth summations in equation (Al). The resulting expan- 
sions were made relative to the nominal conditions selected in the previous 
interpolation. For instance, for the a variable, the regression equations for 
the unknown C, H1 . coefficients were as follows: 


P2 


U, 0, M 0 , j) 


l J 

max J max 


y v v 

j=0 


ijk 


i=0 


k= 1 


(a , 0, M )w* 

Vo °! \ 


(J)dj (J)o k (J) 

bp 


(A3) 


where 

J = 1 • - * n(o, m q ^ 


The quantity, n(o, M 0 ), is the number of data points used to determine 
P 2 (a, 0, M 0 \ as a function of Wj^, 6 sp , and a. 

The function, P 2 (a, °> M o* J ) was computed from the following equation. 

P 2 (a» 0> M 0 , j) = Pi(a w , 0, M q , j) - Pi ^ , 0, M q , (A4) 


The functions, Pi^, 0, M 0 , j) and , 0, M 0 , j), were determined through 

use of the first summation given in equation (Al) and the appropriate Cjj 
coefficients. For instance, the equation for Pl( a w 0 > 0, M 0 , J) is 


max -*max 


Pi 


111CLW IHOA v 

(v °- M o- J ) = EE c ij (v °’ V 


(J)6j(J) 


(A5) 


1=0 J =0 


A similar procedure was used to obtain the expansions for the 3 and M 
variables in equation (Al) . 
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APPENDIX B 


SIMULATION PROGRAM LISTING 


A listing of the computer program for the representation of the YF12 
aircraft motions and propulsion system dynamics is provided in this section. 

A list of the principal fortran variables is given in Table 2. Listings of 
other subroutines used in the execution of the simulation such as an executive 
monitoring program, which provides overall control, and numerical integration 
programs are not included. 

The subprograms shown in the listing are used in the following sequence. 

In subroutine SETUP, the necessary initial quantities are calculated for use 
in subsequent portions of the simulation. Each integration step is controlled 
by subroutine LOOP1. This subroutine calls, in turn, the subroutines for the 
various parts of the simulation - INLET, ENGINE, CONTRL (SAS) , and AIRFRM. 

In each of these subroutines, the necessary constants for the dynamic equations 
are calculated and calls are made to the appropriate integration subroutines 
when needed. Brief descriptions of several additional subroutines which are 
called in the listing are as follows. 

ZXFORM - This subroutine calculates the coefficients to be used in the 
z-transform integration as a function of transfer function coefficients and 
integration interval. 

XFERDA - This subroutine is used in scaling, biasing, and time multi- 
plexing signals to be output to analog recorders. 

MGRATE - This subroutine is the initialization portion of the equation of 
motion integration method. The number of variables to be integrated, order of 
integration, and step size are established. 

GRATON - This subroutine computes and stores results for an integration 
interval for the equations-of-motion integration method. 

ARDC62 - This subroutine is a tabulation of density and speed of sound 
versus altitude for the 1962 ARDC atmosphere. 

VALLMT - This function supplies limits for an indicated variable. 

The computer listing for the aircraft equations is as follows. 
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1 c title setup 

2 subroutine SETUP 
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c 

23 *GIP»P3I 
AsiyP»ANY 

C 

C AlLEPSN ANCLE 

c 

DAT *.^*<PSASP+PSAS> 

° S ASr**PS A S 
C 

if U M SDf> 31# 1*90# 32 

3 1 ! r UCOA.FO.O) S* T? 7 2 
C 

FfcLTAA-nATC 

0A3PP *DELTAA-OAD!C«nT 
n A|P «PSAC»PAADlC*nT 

98 78 33 

c 

32 DELTA4>Z821(2)*DA0P*Zft*?(2M!li)PP«7!2]f ?)*3Af**T??'2>«D4lF 

C 

PASPPOABP 
hair -DAi 
c 

33 DA*P«DELTAA 

c 

C ELEVATOR ANCLE 

c 

nET "*5* 1 05 ASP+OSAS ) 

OSASP-ORAS 

c 

TF UMODH 41# 100# 4? 

41 IF UCDE.EO.O) 08 TS *2 
C 

oeltae*dfjc 

DESPP -DELTAE-D£DIC*3T 
GB T& 43 
C 

42 DFLTAE-2821 (3)*DE6p4?822(3)«3f8P«*^7!?n3)*nFI*7!22t3UDFlP 
r 

DEBPP-DE& P 
C 

43 Df SPwCELTAF 
OE ! P«DE 1 

C 

DETL»DFLTAE*DF7 

C 

C RUDDER ANGLE 

c 

DR I *.5MRSASP*R$AS> 

RSASP-PSAS 

r : 

IF UMSDEl 51*1 00#52 
51 IF UCDR.r0*0) CS T& 52 
C 

DELTAR-DRTC 

drspp «oeltar-droic*dt 
DRIP «RSAS-PSADIC#DT 
06 TA 57 
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156 

m 

15$ 

159 

16C 

161 

162 

163 

164 

165 

166 


C 

52 DeLTAR-Z62U4)«DR0P+7622(4 j«DPf^PP47I21 U)*Dp!«9>7!22(4 UDRIP 

C 

DftflPP“DReP 
OR I P aORf 
C 

53 DR©PotlELT AR 

C 

100 RETURN 
C 

FNO 
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1 

C 


title 



AlRFPM 


2 

3 

c 

SU&POUTUE airfrm 





4 


common /BLK1 

/ i none 

s\*T 

. ITIMES 

.Dy 

.NDT 

5 


1. 

time 

.02* 

.n 




6 

c 








7 


COMMON /PLK2 

/ UB 

. V3 

.W* 

.P? 


6 


1* 

*9 

. PM IR 

♦TMETP 

.PS1P 

. PV 1 1 

.THET 

9 


2# 

PSI 

.DM 

.BFT4 

.ALFA 

.BETA* 

.ALFA* 

10 


3# 

alt 

.ANY 

.AMZ 

*ClP 

.ClCA 

.ClDP 

11 


4 « 

CLDfl p 

,CL^S© 

,CL P R 

#CL p P 

. CKB 

. CNDA 

12 


5# 

Cm Dr 

.CNDBp 

.0*0$* 

.CN*R 

.CN*B 

.peltaa 

13 


6 • 

Deltae 

jDflTAR 

.927 

.$77 


. P T?7 

14 


7 . 

DSPf2> 

#ddp(o> 

> P T2*<2> 

. p LM?i 

# A J C J ,? J 

# PPMM<7) 

15 


9 * 

T|2> 

,NEC1<2> 





16 

c 







17 


COMMON /QLK3 

/ CDZ 

.C!U! 

.C?a? 

*CL 7 

.CLA1 

18 


1# 

CMZ 

. CM A i 

.CM*? 

.c*r*e 


*CnDS p 

19 


2. 

CLF^E 

,ClFDmP 

. CLF^ 

# e l f c ? 

,Gr>A-7 

. V7 

20 


3. 

area 

.XhA$$ 

*CHopn 

.SPAN 

. G 

,XXT 

21 


4# 

X7T 

.XlT 

.XMT 

.X^T 

.X? XX 

. YTYY 

22 


5. 

*I7Z 

.71X7 

.cvq 

.CyTp 

.CYF9 

.CLAP 

23 


6 * 

CM A J 

,CMOP 

.C M^r 


#CL«.. 

.CLPP^ 

24 


7, 

CLPSS 

.CL^A S 

• Cl. *>9 5 

.CLDEP$ 

.Cl rsp 5 

.CN63 

25 


8 . 

CMPpS 

.CN*?S 

.CMflA5 

.CmDPS 

.C^T-nos 

.CNDSPS 

26 

c 








27 


COMMON /PLKU/ ALFAZ 

• ALFAZcv 

. r $°' 

.WJEI 

>MtsP7 

28 


1. 

thft? 

#DE2 

.PBIC 

• CBIC 

.pdi: 

. p MHC 

29 


2. 

theic 

.PSIIC 


. r • E T A I r 

.Altaic 

.ALTir 

30 


3* 

I GOSS' ( 2 » 

.?5 B tt<2) 


. !r. r JF 0 i2i 

>dp*!C(2) 

.?^P3i: (2) 

31 


4# 



• 1CAJCM 

.a jicir-i 

. aj* n? > 

. AJCMfci?! 

32 


5# 

ICWPA^l2),WFAniC<2: 


.“rp^ic<r 

1?»5A5 

.P54SIC 

33 


5* 

PSADjc 

jICCSa* 

.CSaSIc 

*lr=s^ 

.«?« AStr 

, p SOir 

34 


7- 

fsaiic 

.RSAIMC 

• ICtlA 

.CaIC 

.curie 

. icde 

35 


9, 

Drtc 

.DEMr 

.ten 3 

#r»ic 

.np^tc 

.xp 

36 


9# 

?r 

J XM? 

» ALTZ 





37 

38 

39 

40 
4 l 

42 

43 

44 

45 

46 

47 
46 
49 
30 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 


C 

c 

!5|MENM*M 9UFFE*t60) 

C 

tUTA NvA*/10/.|NT9RD/2/ 

C 

C 

IF (JMODn 1.100.2 
C 

C I nTFGQ AT I ON PeUTtNE SETUP 

c 

1 OALt M(JPATE(NVAP#DT # aUFFFR - |MT^Pn # T| MF# frs.pp»j 

C 

c 

CALL A9DC62I 4LTZ.SP?$N%PH6) 

OBA&Z »*5*&*0*</7*V7 
C 

9UAN1 «09ARZ*AREA/XmASS 
0UAN2 -AtFAZR*VZ 
C0AN3 -9PAN/(2.*VZ> 

0UAN4 *Q?APZ*0UAN3 
0UAN5 *Ap£A#$P4N/X IxX 
0UAN6 -APEA*SPAN/y 177 
0UAN7 ’»ape:a*c^ord/x?vy 
GUANO *CH0Rtl/f2.*VZl 

GUANO -OBA?7*OUAN8 
CUANICUOUANS/VZ 
QUAN1 1“X!XZ/X!XX 
QUANl2*X J X2/X I ZZ 
OUAN13-1.-OUAN1 1 OMAN 12 
CUA^14«2,«r5PZ 

c 

Al r A2 «ALFAZP*AL r AZP 

CDA7 *CDZ*C0A£«ALF*7^+rnA7*AtrA2 

CLAZ «clz*clai« ALFaZR 

CM AZ »CHZ*CMA1*AIF4ZR*CMA2* AL r A? 

XZ«0UASJi#C0AZ*0UASl«CtAZ»ALFA7p 

ZZ-C-jANl^CLAZ^CUAK’WCnAl^ALFAZp 

CUAMj7-03AR7*CMAZ 

C 

06 TO 3 
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7Q 

79 

60 

M 

63 

63 

04 

65 

66 
67 
66 
S9 

90 

91 

92 

93 

94 

95 

96 

97 
96 
99 

toe 

101 

102 

103 

104 

105 
10 $ 
107 

106 
109 

uo 

m 

112 

113 

114 

115 

116 
117 
US 
U9 
120 
121 
122 

123 

124 

125 

126 
127 
120 

129 

130 

131 

132 

133 

134 

135 

136 

137 
136 

139 

140 

141 

142 

143 

144 

145 

146 

147 
146 
A 49 

150 

151 

152 

153 

154 

155 


C 

c integration 
c 

2 CAUL GRATeN<aUFFER o U*B*bft,V93jVB*WODpWp4>Pe3,PB,ft9D,G0,RQD*Ra 

PH?D 0 PNjpoT*EDaTH£ 7 R # PSiri*PSt&<»ALTn*ALT > 

C 

c 

p Hf *PHTR/D2B 
TH£? °?MF?R/D2R 
PST ° & S f P /D2P 

0 

UB/9600 
B£TAR°V5/VZ 
BETA agETAR/029 
ALFAR«WP/VZ 
ALFA a ALF AR/D2R 
C 
C 

3 SPM I *S | N C PH I P > 

THFTTL aTHFT7«TNET 

OElAuT«ALT-AlTZ 

c 

TALL A*nc62( ALTi>$P , >?N^o"HPJ 
OSAR «,5*BH<*MVZ*U0I 
FG9AR =G3A9«G9APZ 
C 

AWRpR oALFATR-e-ALFAP 

C 

9UAN15.0QAROA PFA/XMASS 
AWRPP2«AWPPR®4HRPP 
C 

CCA^COZ^COAIoAWPPP^COAp® AWRpP? 

CLA“CtpI + Ct Ai » AXRPR 
0HA«*C^ZPCHAl*4WRPR^cMA2*ANRpP2 
C 

c symmetric cuantities 
c 

PBPSnPDP< 1 ) + DPp ( 2 1 
Dft®S«P3P| 1 WD5P<2>-OUANU 
T S =TUMT<2) 

C 

c asvmkftric quantities 

npPAoDBPt n-0SP(2l 
OSPAaDSPf 1 >«DSP(2> 

?A *T <1 > -T{ 2 ) 

C 

C DRAG AND LIFT IN WIND AVES 
C 

drag wOUANtS* (eDA^CQDE^nELTAr+ODDBPADRPS+cDRSPi^S&S) 
Xtlf r T s OUANl5’»tCLA*CLrDEt»DELTAF^CLP! , !&p*nBpi5^CLFn e P»DS®e 
1 ♦QUANe*CLFQ9*a3! 

r 

C ACCELERATIONS IN BODY AXES 
C 

*l«>Xr*XB-*DRAG+7LIFT*AWppR 

A 2 » 0 L*A^MCYB*BETA*CY 3 R*aFL?A**RUAN 3 *CYRB#Re ) 
A3oZz»2B**LRAGi»AWRPR«XUFT 

A4gj0uAN 3* (OSAR® < ClP*BETA*CLAP*AL Fr A*PrTA*CL'>A®DFLTAA + Cl nP*DELTAR 
l ^CLDRP*D?PA^CUDSP d DSPA ) + 01 JAN4* f ClP^RP + CLfR*®* ) ) 

AftraOUAN** f 00 ARM CNB*B£T A 4-CND A* PEL Y A A ♦CNP J?, *DEL T A P ♦•CnD3p« DRP A 
1 *CNQSP<>J]$PA J*QUAN4MCNP«*PP+CNP3*P^> ) 

GUAN16*A6*XNT*TA 

OUANl 0 — v 7 eRay>G*SPMT+OrjAr 42 *p^ 

C 

UBD a AJ~QU4N2*QB-G*TwFTR*XXTM? 

VRD^Ac^CU AN1 6 
K9:i®A3*VZ* 96 + xZT*?* 

C 

A« “(JUAN?* ( OUANIO^CMAD^NRb^^ua^G® C M Q 3*0B-rUANl 7 + ^? AR* f C m A 
1 <’CMP£ftCElY AE+f MDS^o^DPS > > + X*r 

c 

PRD° f A4^CMJAN 1 1 *9tJAN 16 *XLT*Ta 1 /OUAN 1 3 
QBD^AS+VMTaTS 
RBnaQUANl2*PqD»0UANi 6 
C 

c angle batts er change 
0 
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156 

157 
156 
159 
16C 

161 

162 

163 

164 

165 

166 
167 
166 

169 

170 

171 

1 72 


I P«P3 
r 

C 4 L T T T'jTlF ^ATE 

c 

4 itd*vz*theth_ •n<?p-vB* p ^ t p-w^-su an? 

r: 

C Ny AND *12 ** P G 
C 

4 NYM V^-Q^ANl A + 44,#RPDt/G 
C 

AN2*r4PD^vZ-0Pi)/G 

c 

ion »FTU*N 

C 

END 
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TABLE 2.- LIST OF FORTRAN QUANTITIES 


Airframe 

Quantity 

Fortran 

Units 

Description 

u 

UBD 

m/sec 2 

Forward acceleration 

V 

VBD 

m/sec 2 

Right acceleration 

w 

WBD 

m/sec 2 

Vertical acceleration 

u 

UB 

m/sec 2 

Forward velocity 

V 

VB 

m/sec 2 

Lateral velocity 

w 

WB 

m/sec 2 

Vertical velocity 

p 

PBD 

rad/sec 2 

Roll angular acceleration 

q 

QBD 

rad/sec 2 

Pitch angular acceleration 

r 

RBD 

rad/sec 2 

Yaw angular acceleration 

p 

PB 

rad/sec 

Roll angular velocity 

q 

QB 

rad/sec 

Pitch angular velocity 

r 

RB 

rad/ sec 

Yaw angular velocity 

P(o) 

PBIC 

deg/sec 

Roll angular velocity, I.C. 

q(o) 

QBIC 

deg/sec 

Pitch angular velocity, I.C. 

r(o) 

RBIC 

deg/sec 

Yaw angular velocity, I.C. 

o 

* 

PHID 

rad/sec 

Roll rate 

0 

THED 

rad/sec 

Pitch rate 

* 

PSID 

rad/sec 

Yaw rate 

$ 

PHIR 

rad 

Roll angle 

A6 

■fflETR 

rad 

Pitch angle 


PSIR 

rad 

Yaw angle 

<t> 

PHI 

deg 

Roll angle 

A0 

THET 

deg 

Pitch angle 


PS I 

deg 

Yaw angle 


THETZ 

deg 

Pitch angle, reference condition 

e 

THETTL 

deg 

Total pitch angle 

sin (J) 

SPHI 

ND 

Sine of roll angle 

Ko) 

PHIIC 

deg 

Roll angle, I.C. 

A9(o) 

THEIC 

deg 

Pitch angle, I.C. 

iKo) 

PSIIC 

deg 

Yaw angle, I.C. 
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TABLE 2.— LIST OF FORTRAN QUANTITIES - Continued 


Airframe 

Quantity 

Fortran 

Units 

Description 

AM 

DM 

ND 

Mach number increment 

AM(o) 

DMIC 

ND 

Mach number increment, I.C. 

B 

BETA 

deg 

Sideslip angle 

a 

ALFA 

deg 

Angle of attack 

e 

BETAR 

rad 

Sideslip angle 

a 

ALFAR 

rad 

Angle of attack 

3 Co) 

BETAIC 

deg 

Sideslip angle, I.C. 

a(o) 

ALFAIC 

deg 

Angle of attack, I.C. 

“wo 

ALFAZ 

deg 

Angle of attack, reference condition 

“w 0 

ALFAZR 

rad 

Angle of attack, ref. condition 

h 

ALTD 

m/sec 

Altitude rate 

h 

ALT 

m 

Altitude 

h(o) 

ALTIC 

m 

Altitude, I.C. 

tt/ 180 

D2R 

rad/deg 

Conversion factOT, degrees to radians 

6b Ps 

DBPS 

cm 

i 

Symmetric inlet bypass actuator pos. 

6s Ps 

, DSPS 

cm 

Symmetric inlet spike position 

6 t s 

TS 

ND 

Incremental symmetric engine thrust 

6b Pa 

DBPA 

cm 

Antisymmetric inlet bypass actuator 



position 

6s Pa 

DSPA 

cm 

Antisymmetric inlet spike position 

6 t a 

TA 

ND 

Antisymmetric engine thrust 

C D Ca) 

CDA 

ND 

Drag coefficient 

c L C«) 

CLA 

ND 

Lift coefficient 

C m^ 

CMA 

ND 

Pitching moment coefficient 

y°) 

CDAZ 

ND 

Drag coefficient, I.C. 

c L Co) 

CLAZ 

ND 

Lift coefficient, I.C. 

C m^ 

CMAZ 

ND 

Pitching moment coefficient, I.C. 
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TABLE 2.- LIST OF FORTRAN QUANTITIES - Continued 

Airframe 


Quant ity 

Fortran 

Units 

Des cription 

c Df 

CDZ 

ND 

Zero order term in Cp(ot) 

C D a 

CDA1 

1/rad 

First order term in Cp(a) 

c D a 2 

CDA2 

1/rad 2 

Second order term in C^fa) 

C Lf 

CLZ 

ND 

Zero order term in C L (a) 

Cr 

CLA1 

I/rad 

First order term in C^fa) 

Cmf 

CMZ 

ND 

Zero order term in C m (a) 


CMA1 

1/rad 

First order term in C^Ca) 

W 

(MA2 

1/rad 2 

Second order term in C^fa) 

D/m 

DRAG 

m/sec 2 

Drag in stability axes 

L/m 

XLIFT 

m/ sec 2 

Lift in stability axes 

c »«. 

CDDE 

1/deg 

Aerodynamic derivative 

CD «b P 

CDDBP 

1/cm 

Aerodynamic derivative 

C Dfi 

CDDSP 

1/cm 

Aerodynamic derivative 

°sp 




CL «e ! 

CLFDE 

1/deg 

Aerodynamic derivative 

CL<s bp 

CLFDBP 

1/cm 

Aerodynamic derivative 

C L< * 

CLFDSP 

1/cm 

Aerodynamic derivative 

°sp 




% 

CLFQB 

1/rad 

Aerodynamic derivative 

q v 

QBAR 

N/m 2 

Dynamic pressure 


DQBAR 

N/m 2 

= 

qv 0 

QBARZ 

N/m 2 

Dynamic pressure, reference cond. 

V 

VZ 

m / sec 

Initial forward velocity 

S 

AREA 

m 2 

Wing area 

m 

XjMASS 

kg 

Mass of vehicle 

c 

CHORD 

m 

Mean aerodynamic chord 


i 
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TABLE 2- LIST OF FORTRAN QUANTITIES - Continued 




Airframe 

Quantity 

Fortran 

Units 

Des cription 

b 

SPAN 

m 

Wing span 

P 

RHO 

kg/m 3 

Air density 

g 

G 

m/sec 2 

Acceleration due to gravity 

X 

o 

XZ 

m/sec 2 

Forward acceleration calculated for 
ref. cond. 

Zo 

ZZ 

m/sec 2 

Downward acceleration calculated for 
ref. cond. 

X b 

XB 

m/ sec 2 

Biased quantity for forward 
acceleration 

Zb 

ZB 

m/ sec 2 

Biased quantity for downward 
acceleration 

M b 

XMB 

rad/sec 2 

Biased quantity for pitching angular 
acceleration 

X6 t 

XXT 

m/sec 2 

Forward acceleration due to thrust 

Z ^t 

XZT 

m/sec 2 

Downward accel. due to thrust 

L *t 

XLT 

rad/ sec 2 

Rolling acceleration due to thrust 

M6 t 

XMT 

rad/sec 2 

Pitching acceleration due to thrust 


XNT 

rad/sec 2 

Yawing acceleration due to thrust 

*xx 

XI XX 

kg m 2 

Rolling moment of inertia 

1 yy 

XIYY 

kg m 2 

Pitching moment of inertia 

*ZZ 

XIZZ 

kg m 2 

Yawing moment of inertia 

Ixz 

XI XZ 

kg m 2 

Product of inertia 


CYB 

1/deg 

Aerodynamic derivative 

c y& T 

CYDR 

i 

1/deg 

Aerodynamic derivative 

C Yr 

CYRB 

1/rad 

Aerodynamic derivative 

% 

CLB 

1/deg 

Aerodynamic derivative 


CLAB 

1/deg 2 

Aerodynamic derivative 


1 

CLDA j 

1 

1/deg 

Aerodynamic derivative 
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TABLE 2.- LIST OF FORTRAN QUANTITIES - Continued 




Airframe 

Quantity 

Fortran 

Units 

Description 

C *fi r 

CL DR 

1/deg 

Aerodynamic derivative 

C **bp 

CLDBP 

1/cm 

Aerodynamic derivative 

% 

CLDSP 

1/cm 

Aerodynamic derivative 

Ca 

CLRB 

1/rad 

Aerodynamic derivative 

C* 

CLPB 

1/rad 

Aerodynamic derivative 

Cm a 

CM AD 

1/rad 

Aerodynamic derivative 

c m q 

CMQB 

1/rad 

Aerodynamic derivative 

Cm <5 e 

CMDE 

1/deg 

Aerodynamic derivative 


CMDBP 

1/cm 

Aerodynamic derivative 

S 

CNB 

1/deg 

Aerodynamic derivative 

C n 6a 

CNDA 

1/deg 

Aerodynamic derivative 

Cn 5 r 

CNDR 

1/deg 

Aerodynamic derivative 

^bp 

CNDBP 

1/cm 

Aerodynamic derivative 

S P 

CNDSP 

1/cm 

Aerodynamic derivative 

Cn r 

CNRB 

1/rad 

Aerodynamic derivative 

Cn 

-n p 

CNPB 

1/rad 

Aerodynamic derivative 

% 

CLBS 

1/deg 

Aerodynamic derivative 

c s 

% 

CLPBS 

1/rad 

Aerodynamic derivative 

cl r 

CLRBS 

1/rad 

Aerodynamic derivative 


CLDAS 

1/deg 

Aerodynamic derivative 

r s 

C A<S r 

CLDRS 

1/deg 

Aerodynamic derivative 

C L 

CLDBPS 

1/cm 

! 

Aerodynamic derivative 
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TABLE 2.- LIST OF FORTRAN QUANTITIES - Continued 


Airframe 

Quantity 

Fortran 

Units 

Description 

c ‘*sp 

CLDSPS 

1/cm 

Aerodynamic derivative 


CNBS 

1/deg 

Aerodynamic derivative 


CNPBS 

1/rad 

Aerodynamic derivative 

4 r 

CNRBS 

1/rad 

Aerodynamic derivative 

Cn 6 a 

CNDAS 

1/deg 

Aerodynamic derivative 


CNDRS 

1/deg 

Aerodynamic derivative 

c s 

n *bp 

CNDBPS 

1/cm 

Aerodynamic derivative 

S 

^ n 6 

°sp 

CNDSPS 

1/cm 

Aerodynamic derivative 

cos a Wo 

CALFAZ 

ND 


sin ot Wo 

SALFAZ 

ND 


cos 2 c in, 

0 

CALF Z 2 

ND 


sin 2 a Wo 

SALFZ2 

ND 


n z 

ANZ 

g 

Normal acceleration 


DT 

sec 

Frame time 

SAS Control 

Psas 

PSAS 

deg 

Roll SAS 

<lsas 

QSAS 

deg 

Pitch SAS 

r sas 

RSAS 

deg 

Yaw SAS 


DELTAA 

deg 

Aileron angle 

A6 e 

DELTAE 

deg 

Elevator angle 

6 r 

DELTAR 

deg 

Rudder angle 

«e 

DETL 

deg 

Total elevator angle 

k PP 

XKPP 

sec 

SAS gain 

Psas C“M 

PSASU 

deg 

Psas upper limit 

Psas ^ 

PSASL 

deg 

Psas lower li m it 
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TABLE 2.- LIST OF FORTRAN QUANTITIES - Continued 




SAS Control 

Quantity 

Fortran 

Units 

■ 

Description 

k qq 

XKQQ 

sec 

SAS gain 

k q2 

XKQ2 

sec 

SAS gain 

<tsas 

QSASU 

deg 

Total q sas upper limit 

Qsas 

QSASL 

deg 

Total q sas lower limit 

^sasj, 

QBIU 

deg 

Lagged qsas upP eT limit 

^sas£ 

QBILL 

deg 

Lagged q sas lower limit 

k rr 

XKRR 

sec 

SAS gain 

k m 

XKRN 

deg/g 

SAS gain 

n yn 

ANY 

g 

Lateral acceleration at nose 

r sas ( U ' L '> 

RSASU 

deg 

r sas u PP er limit 

r sas (L.L.) 

RSASL 

deg 

r sas l° wer limit 

Tpi 

TCI 

sec 

Time constant 

^r2 

TC2 

sec 

Time constant 

T q2 

TCN1 

sec 

Time constant 

T ql 

TCD1 

sec 

Time constant 

Tqs 

TCN2 

sec 

Time constant 

T 

<13 

TCD21 

sec 

Time constant 


TCD22 

sec 

Time constant 

T« 

TCN3 

sec 

Time constant 

P 

PBDEG 

deg/ sec 

Roll angular velocity 

q 

QBDEG 

deg/sec 

Pitch angular velocity 

r 

RBDEG 

deg/sec 

Yaw angular velocity 

w n 

n a 

WNA 

rad/sec 

Undamped natural frequency for 
aileron servo 

Pa 

Ik 

ND 

Damping ratio of aileron servo 

«n e 

WNE 

rad/ sec 

Undamped natural frequency for 
e lev on servo 

Pe 

ZE 

ND 

Damping ratio of elevon servo 

Wn r 

WNR 

rad/sec 

Undamped natural frequency for 
rudder servo 

p r 

ZR 

ND 

Damping ratio of rudder servo 
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TABLE 2.- LIST OF FORTRAN QUANTITIES - Continued 


SAS Control 

Quantity 

Fortran 

Units 

Description 

Psas(°) 

PSAS I C 

deg/sec 

!• f° r Psas 

^sas (o) 

QSASIC 

deg/sec 

I.C. for q sas 

« a Co) 

DAIC 

deg 

I.C. for aileron 

M°) 

DBIC 

deg 

I.C. for elevon 

« r (o) 

DRIC 

deg 

I.C. for rudder 

M°) 

DEZ 

deg 

Reference condition for elevon 

Psas C°) 

PS AD I C 

deg/sec 2 

I.C. for p sas rate 

*sas(°) 

PSADIC 

deg/sec 2 

I. C. for r sas rate 

6 a C°) 

DADIC 

deg/sec 

I.C. for aileron rate 

M°) 

DEDIC 

deg/sec 

I.C. for elevon rate 

s r (o) 

DRDIC 

deg/sec 

I.C. for rudder rate 

r sas 

RSAIIC 

deg 

I.C. for first part of r sas equation 

r sas2^°^ 

RSA2IC 

deg 

I.C. for second part of r sas equation 

Inlet 

AM]^ 

DMT 

ND 

Measured increment of Mach number 

«m 

ALFAT 

deg 

Measured angle of attack 

&m 

BETAT 

deg 

Measured sideslip angle 

6 sp 0 

DSPZ 

cm 

Reference condition for inlet spike 

W-i 

lc o 

WICZ 

lb/sec 

Ref. condition for inlet airflow 

(p S /Pt m )c, 

PSZ 

ND 

Compensation for signal pressure at 
reference condition 

(pt /Pt 0 ) 0 

PT2Z 

ND 

Compensation for pressure recovery at 
reference condition 


PSC 

ND 

Command signal for bypass control 

Ssp 

DSP 

cm 

Inlet spike position 

Ps/Pt 0 

PS 

ND 

Bypass signal pressure 

A$bp 

DBP 

cm 

Bypass actuator position 

5 bp 

DBPR 

cm 

Total bypass actuator position 
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TABLE 2.- LIST OF FORTRAN QUANTITIES - Continued 


Inlet 

Quantity 

Fortran 

Units 

Description 

5 bp 0 

DBPZ 

cm 

Ref. condition for bypass actuator 

Aw 4 
A c 

DWIC 

kg/ sec 

Inlet airflow 

Pt2^Pt 0 

^(pt 2 /pt 0 ) 

PT2 

ND 

Inlet pressure recovery 

PT2S 

ND 

Incremental pressure recovery 

f( 8 m) 

G1 

ND 

Spike 8 bypass command function 


G2 

ND 

Bypass command function 

f bp( n %) 

G3 

ND 

Bypass command function 

f b( w i c »*sp»°) 

G4 

ND 

Signal pressure function 

f b( w i c > M ) 

G5 

ND 

Signal pressure function 

f sp( n zm) 

G 6 

in. 

Spike command function 

^ sp (“m) 

G7 

in. 

Spike command function 

^b( w i c >^sp»^) 

G 8 

ND 

Signal pressure function 

fb(wi cJ 6 S p) 

G9 

ND 

Signal pressure function 

6 spd 

G10 

in. 

Unstart boundary function 

f Pa( S sPd)’ fc a^> 

Gil 

kg/sec 

Unstart boundary function 

f =aW 

G12 

kg/sec 

Unstart boundary function 

f c a ( 4 sp.M)*ki 15 4M 

G13 

kg/ sec 

Unstart boundary function 

^a( w i c >^spi a ) 

G14 

ND 

Pressure recovery function 

f a( w ic > 5 sp) 

G15 

ND 

Pressure recovery function 

f a( w i c > 5 sp^) 

G16 

ND 

Pressure recovery function 

f a( w i c sM ) 

G17 

ND 

Pressure recovery function 

f d( w i c ) 

G18 

cm-sec/kg 

Shockwave position function 

fl GO 

G19 

ND 

Incremental press, recovery function 

f c s («) 

G20 

in. 

Unstart boundary function 

T il 

Til 

sec 

Time constant 

T i 2 

TI2 

sec 

Time constant 

T i4 

TI4 

sec 

Time constant 

ki 5 

XKI5 

kg/ cm- sec 

Bypass loop gain 

k i 6 

XKI 6 

cm 

Bypass loop gain 
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TABLE 2.- LIST OF FORTRAN QUANTITIES - Continued 


Inlet 

Quantity 

Fortran 

Units 

Description 

k i 10 

XK110 

ND 

Bypass command gain 

k il2 

XKI12 

cm 

Spike command gain 

k il5 

XKI15 

kg/ sec 

Unstart function gain 

ki 
1 16 

XKI16 

kg/cm-sec 

Unstart function gain 

ki 'l9 

XKI19 

cm 

Spike gain 

ki 2Q 

XKI20 

kg/ sec 

Engine airflow gain 

k, 

lo o 

XKI23 

ND 

Incremental pressure recovery gain 

w n 

“so 

WNSP 

rad/sec 

Undamped natural freq. of spike servo 

p sp 

ZSP 

ND 

Damping ratio of spike servo 

t M 

TIDM 

sec 

Time constant 

T 3 

TIDB 

sec 

Time constant 

T a 

TIDA 

sec 

Time constant 

f bD ( a w n >°) 

G2Z 

ND 

Ref. cond, for bypass command function 

%(“w 0 ) 

G7Z 

cm 

Ref. cond. for spike command function 


ICDSP 

ND 

I.C. mode control for z-transform 




integration 


ICDBP 

ND 

I.C. mode control for z-transfoxm 




integration 

5 Sp^ 0 ^ 

DSPIC 

cm 

I.C. for spike position 

6 sp (°) 

DSPDIC 

cm/ sec 

I.C. for spike position rate 

6 bp (0) 

DBPIC 

cm 

I.C. for bypass position 

6bpC°) 

DBPDIC 

cm/ sec 

I.C. for bypass position rate 


Engine 

PLA 

PLA 

deg 

Power level angle 

A-j 

AJC 

cm 2 

Command signal for exhaust nozzle 

3c 



area 

N 

RPMM 

% 

Command incremental engine rotor 




speed 

N m 

RPM 

% 

Incremental engine rotor speed 

Pt 4 /Pt 0 

PT4 

ND 

Primary burner pressure ratio 
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TABLE 2.- LIST OF FORTRAN QUANTITIES - Continued 


Engine 

Quantity 

Fortran 

Units 

Description 

w fpb 

WFPB 

% 

Incremental primary burner fuel flow 
rate 

w fab 

WFAB 

% 

Incremental afterburner fuel flow rate 

A j 

AJ 

cm 2 

Incremental exhaust nozzle area 

fi t 

T 

% 

Incremental thrust 

e c 

WEC 

kg/ sec 

Incremental airflow demanded by engine 

l T„ 

TE1 

sec 

Time constant 

^e3 

TE3 

sec 

Time constant 

T e4 

TE4 

sec 

Time constant 

^edl 

TED1 

sec 

Time constant 

T en 2 

TEN 2 

sec 

Time constant 

T ed2 

TED2 

sec 

Time constant 

K e3 

XKE3 

cm 2 /%N 

Gain 

K e4 

XKE4 

(kg/sec) /%N 

Gain 

K e ? 

XKE7 

%6^/cm 2 

Gain 

K e8 

XKE8 

ND 

Gain 

K e9 

XKE9 

ND 

Gain 

K el o 

XKE10 

%6 t /%Wfp b 

Gain 

Ken 

XKE11 

1/cm 2 

Gain 

K ei2 

XKE12 

ND 

Gain 

K el 3 

XKE13 

1/%N 

Gain 

K ei4 

XKE14 

%Wf p b/%N 

Gain 

K eis 

XKE15 

%w fpb 

Gain 

K el6 

XKE16 

u t /u fab 

Gain 

K el7 

XKE17 

%w fab /deg 

Gain 

Kel 8 

XKE18 

ND 

Gain 

K el9 

XKE19 

%w fab 

Gain 

K e 20 

XKE20 


Gain 


ICRPM 

ND 

I.C. mode control for z-transform 




integration 

N(0) 

RPMIC 

%N 

I.C. for engine rotor speed 
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TABLE 2 .- LIST OF FORTRAN QUANTITIES - Concluded 


Engine 

Quantity 

Fortran 

Units 

Description 


ICAJ 

ND 

I.C. mode control for z-transform 
integration 

AjCO) 

AJIC 

cm 2 

I.C. for exhaust nozzle area 

Aj (0) 

AJDIC 

cm 2 / sec 

I.C. for exhaust nozzle area rate 

AJCDIC 

cm 2 /sec 

I.C. for command exh. noz. area rate 


ICWFPB 

ND 

I.C. mode control for z-transform 
integration 

w fpb c°) 

WFPB I C 

%w fpb 

I.C. for primary burner fuel flow 

ICWFAB 

ND 

I.C. mode control for z-transform 
integration 


WFABIC 

Wab 

I.C. for afterburner fuel flow 
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Figure !• — Three-view drawing of the aircraft. 
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